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A cavity optomechanical magnetometer is demonstrated. The magnetic field induced expansion of a magne- 
tostrictive material is resonantly transduced onto the physical structure of a highly compliant optical microres- 
onator, and read-out optically with ultra-high sensitivity A peak magnetic field sensitivity of 400 nT Hz~^/^ 
is achieved, with theoretical modeling predicting the possibility of sensitivities below 1 pT Hz~^/^. This chip- 
based magnetometer combines high- sensitivity and large dynamic range with small size and room temperature 
operation. 



Ultra-low field magnetometers are essential components for 
a wide range of practical applications including geology, min- 
eral exploration, archaeology, defence and medicine 111]. The 
field is dominated by superconducting quantum interference 
devices (SQUIDs) operating at cryogenic temperatures 10]. 
Magnetometers capable of room temperature operation of- 
fer significant advantages both in terms of operational costs 
and range of applications. The state-of-the-art are magne- 
tostrictive rnagnetometers with sensitivities in the range of 
fT Hz"^/^!!,!^, and atomic magnetometers which achieve im- 
pressive sensitivities as low as 160 aT Hz"^/^|[5|] but with lim- 
ited dynamic range due to the nonlinear Zeeman eflTectlH, 0] . 
Recently, significant eflTort has been made to miniaturize 
room temperature magnetometers. However both atomic 
and magnetostrictive magnetometers remain generally lim- 
ited to millimeter or centimeter size scales. Smaller mi- 
croscale magnetometers have many potential applications in 
biology, medicine, and condensed matter physics ITI [SD . A par- 
ticularly important application is magnetic resonance imag- 
ing, where by placing the magnetometer in close proxim- 
ity to the sample both sensitivity and resolution may be 
enhanced 13, potentially enabling detection of nuclear spin 
noise imaging of neural networks JtI], and advances in 
areas of medicine such as magneto-cardiography|[ll 0] and 
magneto-encephalographylflll] . 

In the past few years, rapid progress has been achieved on 
NV center based magnetometers. They combine sensitivities 
as low as 4 nT Hz"^/^ with room temperature operation, op- 
tical readout and nanoscale size lll2|] and are predicted theo- 
retically to reach the fT Hz"^/^ range |[l3l] . This has allowed 
three-dimensional magnetic field imaging at the micro scale 
using ensembles of NV-centers and magnetic resonance 
and field imaging [fisl] at the nanoscale using single NV 
centers. In spite of these extraordinary achievements appli- 
cations are hampered by fabrication issues and the intricacy 
of the read-out schemes ifish . Furthermore miniaturization is 
limitied by the bulky read-out optics, the magnetic field coils 
for state preparation and the microwave excitation device! 7]. 

In this letter we present the concept of a cavity optomechan- 
ical field sensor which combines room temperature operation 



and high sensitivity with large dynamic range and small size. 
The sensor leverages results from the emergent field of cavity 
optomechanics where ultra- sensitive force and position sens- 
ing has been demonstrated[16]. A cavity optomechanical sys- 
tem (COMS) is used to sense magnetic field induced deforma- 
tions of a magnetostrictive material, which are detected with 
an all-in-fiber optical system suitable for the telecom wave- 
length range. The presence of mechanical and optical res- 
onances greatly enhances both the response to the magnetic 
field and the measurement sensitivity. 

Three implementations using diflFerent types of COMS are 
investigated, microscale Fabry-Perot resonators IitI] . optome- 
chanical zipper cavities ifisl] . and toroidal whispering gallery 
mode (TWGM) resonators 11 191] . all of which are approxi- 
mately 50 yum in size. Theoretical modelling predicts ul- 
timate Brownian noise limited sensitivities at the level of 
200 pT Hz-i/2^ 10 pT Hz-i/2^ and 700 fT Hz-^/^ respectively 
for each architecture. We experimentally demonstrate the con- 
cept with a TWGM resonator fabricated on a silicon wafer, 
achieving a peak sensitivity of 400 nT Hz"^/^. The possibil- 
ity to integrate TWGM resonators in a dense two-dimensional 
array with waveguides for optical coupling |,2Qi] provides the 
potential for highly-integrated imaging magnetometers 112 1 ] . 



A COMS based field sensor consists of a field- sensitive ma- 
terial, i.e. a magnetostrictive material coupled to the COMSs 
mechanical oscillator. On application of a modulated field 
the field- sensitive material generates an oscillating mechan- 
ical stress field within the COMS structure exciting its me- 
chanical eigenmodes. The vibrations modulate the length of 
the COMSs optical resonator. By coupling light from a laser 
to the cavity at or close to the wavelength of a suitable op- 
tical resonance, the vibrations are imprinted on the transmis- 
sion spectrum and thus can be detected all-optically. To es- 
timate the field sensitivity the COMS based field sensor is 
modelled by a forced harmonic oscillator with a single me- 
chanical eigenmode with eigenfrequency cOm, eflTective mass 
m, and quality factor Q. In this case the displacement of each 
volume element of the COMS u(f, t) = x(t)u(r) is given by 
the product of the scalar displacement amplitude x(t), and the 
spatial mode shape function w(r)||22|]. 
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FIG. 1: (color online). Maximum B -field sensitivity of three COMS 
based magnetometers, a) Cross-section of a TWGM resonator con- 
sisting of a Terfenol-D disc surrounded by an optical resonator (light 
blue), b) Cross-section of a Fabry-Perot resonator with Terfenol- 
D actuators that excite vibrational modes in a suspended micro- 
mirror El. c) Top -view of an optomechanical zipper cavity attached 
to two Terfenol-D actuators fisll . Red arrows represent movement of 
the Terfenol-D (in gray), black arrows mechanical vibrations. 



In the frequency domain the motion of a forced har- 
monic oscillator driven at frequency co is given by x(<jj) 

;\r(^^)[^sig(^^) + ^therm(^^)], whcrC x((^) = [^(^^ 

iojfnOJ / Q)]~^ is the mechanical susceptibility, i^therm(^^) the 
Brownian noise force with the spectral density (l^therm(^^)P) = 
2kTmojm/O i23\\ , and /^sig(^) = ^sig^(^ - ^sig) an effective 



harmonic driving force at the signal frequency cosig- This 
force is generated by the time-dependent mechanical stress 
tensor T(f)e^^''^^ induced by the external magnetic field, which 
produces the body force density /ig(r) = V • T(r)||2l. The 
strength of the effective scalar driving force Fsig experienced 
by the mechanical mode depends on the overlap with the spa- 
tial mode shape function [24, 25] and is given by 



fsigir) • u(f)dV. 



(1) 



The displacement amplitude x(oj) shifts the resonance fre- 
quency of the optical resonator by SD. =_2x(aj), where 
g is the opto-mechanical coupling constant 11221] . In an ex- 



periment this frequency shift is imprinted on the transmission 
spectrum of the optical cavity and measured with a spectrum 
analyzer. Including the spectral density of measurement noise 
S^^2^(oj) due to electronic noise, laser noise and thermorefrac- 
tive noise, the resulting resonance frequency shift spectrum is 
^sa^i^) = ^^(U(^)P) + (^)- Expanding this expression 
we find 



IkTmcOn 



Q 



where the signal field is assumed to be single frequency at 
frequency co^ig. The minimum detectable force F"^^ is ob- 
tained by integrating the signal and noise contributions over 
the bandwidth of the measuring system Ao; and determining 
the force at which the signal-to-noise ratio is unity with the 
result 



2kTmto„ 



cmeas/, ,\ 



Q 



1/2 



Assuming a homogeneous magnetic field oriented in the 
x-direction and a suitably oriented magnetos trie tive medium 
with only a single non-zero magnetos trie tive coefficient dmag 
which causes the material to stretch in the same direction, 
the magnetostrictive induced stress tensor T has only a sin- 
gle component given by Txx = c^magBx- Eq. ([T]) then yields the 
effective driving force that excites the mechanical eigenmode 
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Ux(r)dV = BxC^, 



(2) 



where Ux(r) is the x-component of u(r) and Cact characterizes 
how well the magnetic field is converted into an applied force 
on the oscillator and is referred to here as the magnetic actu- 
ation constant. The minimum detectable magnetic field can 
then finally be expressed as 
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(3) 



The largest reduction is achieved at the resonance frequency 
ojm, where \x(oj)\ is maximized. As expected, the isolation 
from the surrounding heat bath and the enhanced mechanical 
motion achieved by coupling the magnetostrictive material to 
an oscillator enhances the sensitivity. In the thermal noise 
dominated limit, the sensitivity scales as m^^^/CactG^^^, while 
in the measurement noise limited regime it scales as mjc^aQg- 
To estimate the achievable field sensitivity we consider 
three realistic scenarios based on different COMS, a TWGM 
resonator, an optomechanical zipper cavity ifisl] . and a minia- 
ture Fabry-Perot resonator with suspended micro-mirror 
as depicted in Fig. [T^-c. In each case, the magnetostriction 
is provided by the rare earth alloy Terfenol-D which has a 
relatively large magnetostrictive coefficient of a^mag = 5 x 



10^ N T-^m-2 at room temperature 1261 l27l] .The TWGM res- 
onator based scenario is of most relevance to our experiments 
and consists of a central Terfenol-D cylinder of 24yum radius 
and 5yum height surrounded by a cylindrically shaped silica 
ring with 7? = 30 [im outer radius which constitutes the opti- 
cal resonator. We consider a radial breathing mode with typi- 
cal optomechanical coupling constant and mechanical quality 
factor of ^ = Qo/^ = 6 x 10^^ m-^s-^llll and Q = 1000 
respectively, where the displacement amplitude x is defined 
to correspond to the radial displacement at the outer surface 
of the ring. A finite element simulation (COMSOL) yielded 
the eigenfrequency cOm = 27i x 5.6 MHz and effective mass 
m = 35 ng[22], and the magnetic actuation constant Cact = 
0.19 N/T was obtained by solving Eq. (l2l) |[28l] . A measure- 
ment noise level of Sf^^'ico^^^ = 20 Hz/Hz^/^IH was cho- 
sen at an optical resonance frequency of Qq = 1-8 x 10^^ s"\ 
as is typical for TWGM resonators. Fig. [T^ shows the re- 
sulting magnetic field sensitivity as predicted by Eq. ([3]), with 
a maximum predicted sensitivity of 700 fT Hz"^/^ at reso- 
nance. By comparison, the peak sensitivities predicted for 
miniature Fabry-Perot resonators and optomechanical zipper 
cavities of similar spatial dimensions were 200 pT Hz"^/^ and 
10 pT Hz"^/^, respectively (see Supplemental Information for 
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FIG. 2: (color online). Experimental set-up showing a TWGM res- 
onator with a Terfenol-D sample attached. 



calculations 112 811 ). The Fabry-Perot resonator based system 
was least sensitive, as the magnetostrictive material excites 
the mechanical modes of the suspended mirror indirectly and 
inefficiently by vibrating the entire mirror mount ll28l] . The 
TWGM magnetometer shows the best sensitivity as the mag- 
netostrictive material couples to the vibrating mode over its 
entire radius which results in the largest actuation constant 
Cact of the three systems. 

As a first proof of principal demonstration of a cavity op- 
tomechanical magnetometer, a piece of Terfenol-D (Etrema 
Products Inc.) with a size of roughly 50x15x10 yum^ was af- 
fixed to the top surface of a TWGM resonator using microma- 
nipulators and two-component epoxy. The TWGM resonator 
had major and minor diameters of 60 and 6 yum, respectively, 
and a 10 jim undercut! 19]. Fig. [2] shows the experimental 
set-up. The TWGM resonator was placed between two 20 
mm diameter coils that generated an HE magnetic field previ- 
ously calibrated with a commercial Hall sensor. 980 nm light 
from a widely tunable, external cavity diode laser was passed 
through a fibre polarization controller and evanescently cou- 
pled to the resonator via a tapered optical fibre 1I23I1 . The laser 
was thermally locked to the full width half maximum of an 
optical resonance 1I29I] . Mechanical vibrations of the TWGM 
resonator shift its optical resonances and thus modulate the 
transmitted light that is detected with an InGaAs photodiode. 

The transduction spectrum shown in Fig. [3^ was obtained 
by analyzing the transmitted light with a spectrum analyser. 
It shows multiple characteristic peaks corresponding to ther- 
mally excited mechanical modes. With a finite element solver 
the three vibrational modes around 10 MHz were identified 
as the lowest order flexural mode (centre) and the two lowest 
order crown modes. Application of a 250 yuT magnetic field at 
oj^Qf = 271 X 10.38 MHz resulted in a signal at that frequency 
above the Brownian noise confirming the ability of the sensor 
to detect magnetic fields. The magnetic field sensitivity was 
determined by 

Bf\aj,,f)/Aaj'^^ = \BA/^JSNR(oJr,f) Aoj^bw, (4) 

where SNR is the signal to noise ratio, the applied mag- 
netic field in direction x, and AccJrbw the resolution bandwidth 
of the spectrum analyzer. 

To determine the sensitivity over a wide frequency range. 
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FIG. 3: (color online), (a) Brownian noise spectrum with magnetic 
excitation at 10.38 MHz. (b) System response as a function of ap- 
plied B -field frequency. Inset: magnified system response centered 
around 10.385 MHz. Red curves: Lorentzian fits. 



a network analyzer was used to scan the frequency of the 
magnetic field from 7 to 13 MHz and simultaneously detect 
the systems response N(oj). As expected. Fig. [3]3 shows a 
Lorentzian line shaped peak distribution similar to the Brow- 
nian noise spectrum (Fig. [3^), with an additional peak present 
around 11.5 MHz, due to the enhanced noise rejection of 
the phase sensitive network analyzer. In both cases the mea- 
sured data agrees quite well with a model with five mechanical 
resonances 1I23I] . 

Importantly, the apparent large noise amplitudes visible in 
Fig.[3t) are stationary in time, and in fact corresponds to an en- 
semble of resonant features having Q-factors on the order of 
10^. Two such features are shown in the inset. We believe that 
these resonances are due to ultrasonic waves in the Terfenol- 
D grains that are excited by a magnetostrictive mechanism. 
Similar resonances have been reported for other magnetostric- 
tive materials with comparable g-factorsfl^, EH], though to 
our knowledge, this is the first such observation in Terfenol- 
D. Constructive and destructive interferences between neigh- 
boring resonances cause strong positive and negative varia- 
tions in the spectrum. It is expected that the density of reso- 
nances could be greatly reduced by engineering the geometry 
of the Terfenol-D. For example, the simple geometry used in 
Ref. |[3Q|i meant that only roughly two resonances per MHz 
were observed in agreement with theoretical predictions. 

The sensitivity spectrum B^^((jj) shown in 
Fig. m was obtained using the relation 

(0j)N(aJr.f)/N(aj)S (0Jrcf)Bf^(C0r,f) 



B^^ioj) = 



with the reference 
sensitivity obtained from a measurement similar to that in 
Fig. [3^ using Eq. Due to the Terfenol-D resonances, the 
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FIG. 4: (color online). B -field sensitivity as a function of frequency. 



sensitivity varies between 2yuT Hz"^/^ and 170 nT Hz"^/^, 
depending on whether destructive or constructive interference 
occurs. At resonance we achieve a peak sensitivity of 
400 nT Hz"^/^ . As expected, this value is substantially above 
the 700 fr Hz"^/^ predicted for an ideal system. In our first 
prototype the response of the Terfenol-D was transduced to 
the resonator through an epoxy layer and its geometry was 
not well controlled. As a result the overlap between the 
applied force and the material motion is expected to be quite 
poor. Furthermore, the mechanical modes best coupled were 
crown modes, which themselves do not couple with high 
efiftciency to the optical field|0]. These limitations could be 
overcome by controlled deposition of Terfenol-D via sputter 
coating. Further sensitivity enhancement may be achieved by 
using TWGM resonators with higher compliance such as is 
achieved with spoked structures [3 2]. 

In conclusion, a cavity optomechanical magnetometer is 
reported where the magnetic field induced strain of a mag- 
netostrictive material is coupled to the mechanical vibra- 
tions of an optomechanical resonator. A peak sensitivity 
of 400 nT Hz"^/^ is achieved, with theoretical calculations 
predicting sensitivities in the high fT Hz"^/^ with an opti- 
mized apparatus. This optomechanical magnetometer com- 
bines high- sensitivity, large dynamic range with small size 
and room temperature operation, providing an alternative to 
NV center based magnetometers for ultra- sensitive room tem- 
perature applications, low-field NMR and magnetic resonance 
imaging, and magnetic field mapping with high-density mag- 
netometer arrays. 
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